Higgs production at NNLO in QCD: the VBF channel 
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We present a phenomenological study for the production of the Higgs boson at next-to-next-to-leading order 
(NNLO) in QCD via the vector boson fusion (VBF) process. After a general discussion about the different 
production channels of the Higgs, we show results for hadron colliders like LHC and Tevatron in VBF. The 
theoretical predictions are obtained using the structure function approach. This approximation turns out to be 
more accurate than the precision to which the VBF Higgs production channel can be considered a well defined 
process by itself and the theoretical uncertainty are of the order of 1-2%. The uncertainties due to parton 
distributions are also discussed and are estimated to be at the same level. 



1. INTRODUCTION 

To discover whether the Standard Model Higgs 
boson exists or not is one of the aims with most 
priority at the nowadays hadron colliders. These 
are the proton-anti-proton Tevatron machine at 
Fermilab and the proton-proton LHC machine at 
CERN. To achieve this goal an estimate of the 
expected events and the control over those pro- 
cesses which represent the background noise is re- 
quired. Nevertheless in the case of a discovery of 
the Higgs it would become necessary to also in- 
vestigate precisely its properties. On its turn this 
demands a determination of the cross section in 
the various production channels as precise as pos- 
sible. 

At hadron colliders the Higgs production chan- 
nels which have a large enough cross section to 
be relevant are the gluon-gluon fusion Fig[TlJa), 
the VBF Figdjb) and the associated production 
with W, Z bosons and tt FiglTfc,d). 

Generally the dominant production mechanism 
is represented by the gluon fusion mediated by a 
top quark loop. For this production channel the 
first NNLO corrections in QCD have been com- 
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Figure 1. Higgs boson production channels at 
hadron colliders. 



puted in the context of the effective theory in the 
limit large top quark mass limit |1|2|3|4| . It is also 
known that this approximation works very well 
up to the NNLO. Indeed the impact of the finite 
top mass effects are about a factor of ten smaller 
than the uncertainty due to the scale variation at 
the same order |5|6|7j . As far as the associated 
production of the Higgs with W, Z bosons is con- 
cerned, the NNLO QCD corrections have been 
implemented in [S]. Thanks to these results the 
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theoretical uncertainty is reduced to about 10% 
for the gluon fusion production mechanism and 
to less than a few percent for the associated pro- 
duction with gauge bosons. Very recently NNLO 
corrections in QCD have been included also for 
the VBF production mechanism [9] via the struc- 
ture function approach thus reducing the theo- 
retical uncertainty for this channel from the 5- 
10% of the NLO QCD and electroweak combined 
computations jlOlllj down to 1-2%. This makes 
VBF the theoretically most accurate Higgs boson 
production channel at hadron colliders. 

2. THE VBF SIGNAL 

Among the production mechanisms of the 
Higgs boson that we have mentioned above, the 
VBF production mechanism can be considered by 
itself as a signal where its "background" includes 
also Higgs boson production by mechanisms other 
than VBF (see FigH]). 
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Figure 2. Higgs production cross sections at the 
LHC for the various channels (see Ref.[T2] for 
more details about this plot). 



Here we briefly describe its characteristics. Due 
to the preference of the Higgs to couple with 
heavy particles the VBF at LHC is second in size 
only to the gluon fusion channel which is medi- 
ated by a top quark loop (see Fig 12]). Further- 



more it provides a clean experimental signature: 
it usually consists of two almost back to back hard 
tagging jets (generally with an invariant mass big- 
ger than 600 GeV and a pseudo-rapidity separa- 
tion between them Arj > 4) and the Higgs decay 
is confined in the central rapidity region. Impo- 
sing these constraints to the invariant mass and 
the rapidity of the jets as additional cuts (usually 
called VBF cuts) one reaches an impressive im- 
provement of the signal-to-baekground ratio |13j . 

An heuristic way to qualitatively understand 
the property of the VBF signature is to consider 
angular ordering which arises from the dominant 
soft gluon emissions. One finds that the averaged 
azimuthal emission is confined to a cone limited 
by the angle between the emitter and the specta- 
tor parton. It can easily be shown that because of 
the behavior of the propagators in the i-channel 
in Figlljb) small values of the emission angle are 
favored. Because of this the hadronic activity in 
the VBF production mechanism is concentrated 
in the forward/backward region with respect to 
the incoming colored partons. This is in contrast 
to the gluon fusion channel where the radiation is 
expected and found in the central rapidity region 

Keeping in mind that what we strictly call VBF 
signal is the Higgs boson weak production with a 
color singlet exchange in the t-channel, we want 
now to discuss the possible interference effects 
with other processes. Already at LO the VBF 
process in Fig[TJb) can interfere with the produc- 
tion mechanism in Figlljc) where the associated 
produced gauge boson decays into two quarks. 
However as it is shown in [11] this interference 
effect is at the per mil level. At higher orders, in- 
terference effects occur with the gluon fusion pro- 
duction mechanism in Fig[TJ Also in this case the 
interference effect is small [15|16j and it turns out 
to be well below the percent level. This demon- 
strates that the Higgs boson production via the 
VBF mechanism can be defined a signal by itself 
within an ambiguity better than the 1%. It is this 
ambiguity in defining the Higgs boson production 
via VBF that sets also the target theoretical pre- 
cision for this observable. 



3 



3. THE COMPUTATION 

Here we want to illustrate the structure func- 
tion approach in the VBF production chnnel [TU] 
and show that this approximation remains suffi- 
ciently accurate even at NNLO in QCD. 

The structure function approach consists ba- 
sically in viewing the VBF process as a dou- 
ble deep-inelastic scattering (DIS) attached to 
the colorless pure electroweak vector boson fu- 
sion into a Higgs boson. A qualitative illustra- 
tion of the structure function approach is given 
in Figini According to this approach one can in- 
clude NLO QCD corrections to the VBF process 
employing the standard DIS structure functions 
Fi{x,Q'^); i = 1,2, 3 at NLO [T^. Similarly at the 
NNLO level one has to employ the corresponding 
structure functions jl8|19|20|21] . 




Figure 3. Schematic view of the structure func- 
tion approach. 



The structure function approach represents a 
very accurate approximation because it is based 
on the absence or smallness of the QCD inter- 
ference contributions between the two inclusive 
final states Xi and A2. We now discuss the vari- 
ous contributions up to NNLO which in principle 
violate the structure function approach but which 
can nevertheless be safely neglected. 

At LO there is already a structure function 
violating contribution coming from the interfe- 
rences between identical final state quarks (e.g. 
uu — > Huu) or between processes where either a 
W or & Z can be exchanged (e.g. ud — > Hud). 
Simple kinematical arguments show that such 



contributions are very small and contribute to the 
total cross section well below the percent level 
[35] . These contributions can be easily computed 
and have been included in our results anyway. 

At the NLO level possible contributions viola- 
ting the structure function approach arise when a 
gluon in the i-channel is exchanged between the 
upper and the lower quark line in Figl^Kb). How- 
ever the interference of such one-loop contribu- 
tions with the LO diagram have a vanishing color 
factor due to the generators <° (a = 1, . . . , N'^ — 1) 
of the color group SU (Nc) being traceless. This 
means that apart from the interference effects di- 
scussed at LO the structure function approach re- 
presents an exact approach to the computation. 




Figure 4. Examples of neglected Feynman dia- 
grams at NNLO. 



At NNLO the structure function approach is 
not exact but it can be still considered a very 
good approximation. The types of diagrams that 
violate the structure function approach are shown 
in FigHl The first type represents a double gluon- 
exchange in the t-channel (note that one or both 
of the two emissions could also be real); the se- 
cond type is a representative of the so-called sin- 
gle quark line (SQL) diagrams contributing at 
NNLO; the last two are heavy quarks (top and 
bottom) loop diagrams. The first type of con- 
tributions represents a gauge invariant, infrared 
and ultraviolet finite class of diagrams. Another 
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Table 1 

Total cross section (pb) from the neglected triangle diagram in FiglH In parenthesis are reported also the 
numbers included in the NNLO computation in the structure function approach employing the MSTW 
PDF set US]. Ah the parameters have been taken from PDG 2008 [17]. 



characteristic of this class of diagrams is its color 
suppression by a factor of l/N"^ with respect to 
the contributions included by the structure func- 
tion approach. Furthermore this type of contribu- 
tions are also strongly kinematically suppressed 
[23|24|25] . This is mainly due to the behavior 
of the gluon propagator in the i-channel and/or 
to the small overlapping of the phase space of 
real emissions from the upper quark line and real 
emissions from the lower one. The neglected SQL 
type contributions in Fig|4] do not represent a 
class of infrared safe diagrams. However as shown 
in jl3j their impact is small enough not to pro- 
duce a significant deterioration of the VBF sig- 
nal. Also, these color exchange effects are, by our 
definition, no VBF processes. Finally we take 
into consideration the triangle and the box con- 
tributions in Fig|4l Even if a full computation 
is in progress [35] as a first rough estimation we 
have computed the triangle contribution in Fig[4] 
in the limit of infinite top mass. In Table [1] we 
report some values of the contribution to the to- 
tal cross section from the triangle at 1.96 TeV for 
the Tevatron and at 7 TeV for the LHC. This has 
been checked performing two independent com- 
putations. As we can see from Table [1] its impact 
is very small and can be safely neglected. 

4. RESULTS AT HADRON COLLIDERS 

Here we discuss some numerical results ob- 
tained for the VBF production mechanism at 
hadron colliders up to the NNLO in QCD em- 
ploying the structure function approach. 

For the electroweak parameters we have used 
the central values released by the PDG collabo- 
ration in 2008 [57] while for the parton distribu- 
tions functions we choose the set MSTW [Ml . For 
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Figure 5. Total cross section for the LHC at 7 TeV 
employing the MSTW PDF set [26]. 



the LHC we set the center of mass energy to be 
VS' = 7 TeV for the LHC and ^/S^ = 1.96 TeV for 
the Tevatron. 

In FigOwe plot the total cross section for the 
VBF production mechanism at the LHC at 7 TeV. 
The LO, NLO and NNLO resuhs in QCD are 
shown 0jS 0j function of the Higgs boson mass 
which first of all induce only a rather mildly de- 
pendence on it. The bands represent the theo- 
retical uncertainty of the prediction. They have 
been obtained varying the factorization and the 
renormalization scales in the quite large range 
A^R(F) G [Q/4, 4(5] where Q is the virtuality of 
the vector bosons which "fuse" into the Higgs. 
Clearly other scale choices are possible (e.g. the 
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Figure 6. The same as Figl5]for the Tevatron. 
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Figure 7. The PDF uncertainty of the total cross 
section at NNLO for LHC at 7 TeV for the 68% 
CL MSTW PDF set [26]. For ABKM [29] and 
JR09VF [30] the ration of the central value is 
plotted. 



choice Q — nih the Higgs mass) but the one cho- 
sen for the plot turned out to be the more natural 
choice because it exhibits a better convergence of 
the perturbative expansion. This also shows that 
at NNLO in QCD the theoretical uncertainty is 
reduced to be less than the 2% reaching the same 
level of ambiguity at which the Higgs production 
signal via VBF can be defined phenomenologi- 
cally. 

In Figini we report on numbers for the Teva- 
tron where the center of mass energy is set to 
1.96 TeV. As one can expect the total cross sec- 
tion shows the same behavior upon varying the 
Higgs mass and it is almost an order of magni- 
tude smaller. Again the lower part of Fig [5] shows 
a very good convergence of the perturbative QCD 
expansion. Even if the theoretical uncertainty re- 
mains slightly bigger compared to the LHC the 
relative improvement with respect to the NLO 
prediction is at the same order of percentage. 

Finally we consider also the uncertainties co- 
ming from the parton distributions. To achieve 
this we have employed the MSTW 68% confi- 
dence level PDF sets [55] and compare with other 
NNLO PDF sets, i.e. ABKM gH] and JR09VF 
[3D]. The results in Figs (T] |S] show that an almost 



constant 2% PDF uncertainty can be associated 
to the cross section for both the LHC and the 
Tevatron. In the case of the Tevatron the diffe- 
rence between the MSTW, ABKM and JR sets is 
due to larger uncertainties for the high-x quark 

PDFs [ini- 

To conclude we give the address of the web in- 
terface where our code for the NNLO VBF to- 
tal cross section can be used online [31]. After 
the registration, setting on a dialog window the 
energy and choosing the hadron collider the nu- 
merical answer is received per email. 

5. CONCLUSIONS 

We have described recent progress for the VBF 
Higgs predictions and investigated at which level 
it can be considered a well defined process by it- 
self. Then after showing how well the structure 
function approach works even at the NNLO in 
QCD we have employed it to obtain predictions at 
the hadron colliders. The theoretical uncertainty 
is estimated to be less than 2% which is compa- 
tible with the ambiguity at which this signal can 
be defined. Finally what has been presented here 
is a natural first step towards less inclusive (e.g. 
rapidity distributions) predictions of this process 
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Figure 8. Same as Figl7]for the Tevatron. 



at NNLO in QCD. 
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